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The spatial inhomogcncities are expected to affect nucleation process in an essential way. 
These effects are studied theoretically by considering the case of the depinning of the charge 
density wave as a typical example. The threshold field of the depinning of the one-dimensional 
commensurate charge density wave with one impurity has been examined based on the phase 
Hamiltonian at absolute zero. It is found that the threshold field is lowered by a finite amount 
compared to that in the absence of an impurity. 
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§1. Introduction 

Nucleation is one of the most drastic phenomenon in various fields of physics, chemistry, biology, 
and also in engineering!!!) Especially, the nucleation in condensed matter physics is most interesting 
in the sense that it can be controlled by such parameters as pressure, temperature, electric and 
magnetic fields so on. In general nucleation is defined as a phenomenon where a new phase appears 
locally in space. 

The theoretical analysis of nucleation was given by Langer,& who investigated the problem of 
reversing of the direction of magnetization in a ferromagnetic system. In the process of reversing, 
changes of magnetization is found to be not uniform in space but it is triggered by the appearance 
of the magnetic bubbles called droplets. While this theory offers the fundamental understanding of 
the nucleation in a homogeneous medium, we know that a local inhomogeneity in actual systems 
plays important roles in the nucleation process, which will be studied in this paper. For explicit 
studies we choose the case of the charge density wave (CDW)i> depinning clS cl typical example, 
because the field theory which is described by the phase Hamiltonian has been established $ and 
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the concept of the nucleation can be defined clearly in this case as the appearance of the spatially 
local non-uniform structure of the phase variable which triggers the depinning. 

CDW is one of the characteristic state of quasi-one-dimensional conductors, where translational 
symmetry is spontaneously broken. The gapless sliding mode& that incarnates Goldstone mode 
can be pinned by the external objects which break the translational symmetry. If the external 
object is the underlying lattice whose spatial periodicity is commensurate to CDW's, it is called 
commensurate CDW. In this case CDW behaves as an insulator. However, if electric field is applied, 
CDW starts to move above some critical field called threshold field. To estimate the threshold field 
theoretically, the phase Hamiltonian approach! is useful. With the phase Hamiltonian approach, 
the classical threshold field can be derived as the field at which the potential barrier for the sliding 
disappears. In this paper we investigate classically the depinning processes of one-dimensional 
commensurate CDW with inhomogeneity at absolute zero. 

In §2 we introduce our model based on the phase Hamiltonian. In §3 we review briefly in a 
homogeneous case. Below the uniform depinning field, the nucleation requires a finite excitation 
energy. In §4 we examine the ground state in the presence of an impurity and in §5 we investigate 
the threshold field in this case. The threshold field can be smaller than the uniform depinning field 
and depinning sets at the impurity site. In §6 the potential curve of our model is shown. Our 
conclusion and discussion are given in §7. The effects of three dimensionality and fluctuations, 
quantum and thermal, will be studied separately. 

§2. The Model 

We investigate the one-dimensional commensurate CDW with one impurity located at Xj. The 
Lagrangian is then given by 



dX 



A (j^j -F<f> + g(i- cos(M0)) - V t cos( 2 £: F X + </>)8(X - X { 



(2.1) 



Here, the first term in the integration is the elastic energy, A = (hvp) / (&ir) , with vp being the 
Fermi velocity. The second term is the energy associated with the electric field, F = (eE)/(2ir), 
with E and — e (e > 0) being electric field and the electron charge, respectively. The third term is 
the commensurability energy, g = (| A\ 2 /ep)n e (\A\/W) M ~ 2 , \A\ being the Peierls gap, and ep, n e , 
W are the Fermi energy, the density of electron, and the width of the original band, respectively. 
Here M is the degree of commensurability, M = 7r/(kpa), kp is the Fermi wave number, a is the 
lattice spacing. The last term is the coupling energy of CDW to one impurity. In this Lagrangian, 
there exists a characteristic length, £ = \f (2 A) / (M 2 g) , which is the phase coherence length due 
to the commensurability. We make the Lagrangian dimensionless through scaling by this £; i.e. 
x = C = M 2 g£L, where L is given by 



L = — dx 
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(2.2) 



2 



Here e = F/(gM 2 ), v = V{/ '(g^M 2 ). In this paper, we assume e > 0. Further, we define xi = 
and x = (2nz)/M, where z = Xi/a characterizes the location of an impurity relative to the site 
where the energy gain by commensurability is maximum. The range of x is — vr/M < x ^ tt/M, 
because the Lagrangian has the periodicity of 2-k/M with respect to the phase. 

§3. The Homogeneous Case 

First of all, we examine the case of v = 0, namely the homogeneous case.& 
The stable configuration of the phase is determined by varying Lagrangian; 

-/-e + -^sin(M$=0. (3.1) 

When s = 0, eq. ([O]) is the sine-Gordon equation, and has two kinds of solutions, which are a 
trivial one, <j) = 0, and the kink solution, which is given by 

4 

4>(x) = — arctan [exp(±(x — xo)j = 4>±{x — xo), (3-2) 
and shown in Fig. 1. Here xq is the center of the kink. Among these solutions, the lowest energy 



Fig. 1. The classical solution in the case of e = 0. The solid line is a kink and the dotted line is an anti-kink. 



solution, namely, the ground state solution is </> = 0, and the kink and the anti-kink are excitations 
with finite excitation energy. 
In the case of e 7^ the equation of the classical configuration has also an uniform solution, 

4> = arcsin(eM), (3.3) 
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and the non- uniform solution with local spatial variation, <f>i(x), which is given by 



1 



M 



1 



£ ( — arcsin(eM) - ^ ) + — ( ^/l - ( £ M) 2 - cos(M0, 



(3.4) 



This equation is derived by the integration of eq. (|3.1| ) with the boundary condition, ^(oo) = 
(1/M) arcsin(eM) and 4>'(oo) = 0. The solution, (j)i(x), is shown in Fig. 2. We notice that the 
non-uniform solution, <pi(x), is kink pair solution. 



Fig. 2. The non-uniform solution, (f>i(x), of the classical equation in the case of e = 0.1. 



In the presence of the electric field, the cosine-type potential for the uniform phase is tilted as is 
shown in Fig. 3. The potential barrier disappears at e = 1/M = et, which is the classical depinning 
field of the uniform CDW. 

On the other hand, if the kink pair is excited, CDW will also depin. This is another type of 
depinning process which sets locally in space and considered to be the nucleation process. However, 
below the threshold field, et, of the uniform depinning, the kink pair is always the excitation, which 
requires a finite excitation energy. Therefore the nucleation does not take place in a homogeneous 
case without quantum nor thermal fluctuations. 

§4. Ground State in the Presence of an Impurity 

Next, we consider a case with an impurity, jj ^ 0. The stable configuration of the phase is now 
given by 

- 4>" -e + — sin(M0) + v sin(x + <j>)6(x - Xi) = 0, (4.1) 
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Fig. 3. The potential for the uniform phase at e = 0.1 for a choice of M — 4. 



with the boundary conditions 



(±oo) = — arcsin(eM) , 0'(±oo) = 0. (4.2) 



In the case of e = 0, the solution is given analytically. Depending on the range of x> there are 
two kinds of configurations. If x is in the range of —tt/M < \ < 0, the solution is located in the 
region between and tt/M (Config.l) as shown in Fig. 4. If x is 111 the range of < x < tt/M, 
the solution is located in the region between —tt/M and (Config.2). We take electric field e > 0. 
Under this condition, Config.l has more tendency to depin than Config.2. Namely, in the case of 
—tt/M < x < 0, the threshold field can be smaller than that of uniform dep inning as is disclosed in 
§5; in the case of < x ^ tt/M, however, the threshold field is same as that of uniform depinning. 
Therefore, we consider only the case of —tt/M < x < 0. This choice is justified by the discussion 
in §7. 

The solution of Config.l is obtained in terms of 4>-{x) given in eq. (|3.2|), 



4 

4>c(x) = — arctan [exp (-\x - Xi\ - c )] (4.3) 
= (j)-(\x-Xi\+c ), (4.4) 



where cq is the parameter which is determined by the equation, 

4 1 / 4 \ 

— — — ^- + v sin [x + -tt arctan exp(-co) =0. (4.5) 
M cosn(co) \ M J 

This equation is due to the requirement that the jump of the first derivative of 4> c {x) at should 

match the third term in the left hand side of eq. ((OJ). The solution (j)c{x) is considered to be a 
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Fig. 4. The two kinds of the configurations with respect to x- This figure is for a choice of M = 4, v = 1 and e = 0. 



kink and an anti-kink connected at the impurity site as shown in Fig. 5. Further the parameter cq 
is equivalent to the distance from the impurity site to the center of the kink as indicated in Fig. 5. 

In the case of e ^ 0, the stable solution, (j) s (x), in the presence of an impurity given by eq. Q4.1| ) is 
obtained by the connection at the impurity site of two non-uniform solutions, one of which, 4>i{x), 
has been shown in Fig. 2. An example of the stable solution, (fi s (x), is shown in Fig. 6. 

Though the stable solution, (j) s (x), is obtained only numerically, the information of the phase 
value at the impurity site is obtained. The stable solution, (j) s (x), is given as follows; 



(f) s (x) = (f>i(\x - Xi\ + c) 



(4.6) 



where <f>i(x) is the non-uniform solution of eq. (3.1) and c is a parameter to be determined to satisfy 



the boundary condition given by eq. (4.2). When we consider 4> s (x), it is necessary to consider only 
the range of solid line in Fig. 6. 

Substituting eq. (|4.6|) into eq. (|4.1[), we obtain 



2(f>' l (c)+vsm(x + <t>l( c )) = 0. 



(4.7) 



By use of eq. (|3.4|), eq. (|4.7|) is rewritten as 



e ( arcsin(eM) - 0,(c) J + ^ ( V 1 ~ ( eM ) 2 ~ cos(M# ( c)) 



1 

M 2 



2 1 

+ 2 wsin (x + <M C )) = °- 
(4.8) 

This equation determines the phase value, (f>i(c), (instead of c) at the impurity site of the stable 
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Fig. 5. The solid line is the stable solution with an impurity at i = i, for a choice of M = 4, v = 1, x = —Tf/M 
and e = 0. 



solution when the electric field, e, is given. It means that instead of the boundary condition at 



§5. Threshold Field in the Presence of an Impurity 

The depinning threshold field, e c , is determined as the field at which the stable solution becomes 
unstable^ The instability of the solution is triggered by the onset of a zero eigenvalue of the 
fluctuation mode around the stable solutionis 

The eigenvalue equation of the fluctuations, 5(j)(x), around the stable solution, <j) s {x), is the second 
variational equation of the Lagrangian; i.e. 



where A is the eigenvalue. Hence we consider the case of A = to determine the threshold field, e c . 

First, we consider the fluctuation around the solution, (fti(x), in the case of v = 0. The equation 
of the zero-mode fluctuation, 5(f>i(x), around (j>i(x) is 



infinity, eq. ([4.2|), we obtain a boundary condition at the impurity site. 






(5.2) 



and the boundary condition is 




(5.3) 



Its solution is 




(5.4) 
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Fig. 6. The solid line is the stable solution with an impurity for a choice of M = 4, v — 1, \ = —n/M and e = 0.1. 



Next, we consider the case of v ^ 0. The threshold field, e c , is determined as a particular value 
of the electric field at which zero-mode fluctuations in both sides of the impurity are connected. 
The zero-mode fluctuation around the stable solution in each region must be expressed as 

6<j)(x) oc -{29{x - Xi) - l)(j)' s (x) (5.5) 
= —<// l (\x — x i \+c), (5.6) 

by noting 

<j>' s {x) = {2e{x-x i )-l)4>' l {x), (5.7) 
where 9(x) is the step function. Substituting eq. (5.6) into eq. (|5.1|), we obtain 



2#(c) -vcos( X + HcM(c) = 0. (5.8) 



At the threshold field, e c , eq. (^8) should be satisfied. By use of eqs. ( pH] ) and ( |4,7j ), eq. ( plq ) is 
expressed as follows; 

-e + sin(M^(c))) - ^ 2 sin [2( X + ^i(c))] = 0. (5.9) 

Hence the threshold field, e c , is determined as a value of e at which both eqs. (fh^) and (^^) have 
a common solution, 0/(c), for each fixed values of x and u. 

The solutions of eqs. ( [4.8D and (|5.9[ ) for M = 4 are shown in Fig. 7. The threshold field, e c , 
is normalized by the threshold field, et, of the uniform case. When v is fixed, the normalized 
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Fig. 7. The threshold field in the presence of an impurity, normalized by that of the uniform case, for a choice of 
M = 4. The threshold field for a choice of v = 1 is shown in (a) and for a choice of x — —n/M is shown in (b). 

threshold field goes to 1 as x tends to —n/{2M). It is noted that the threshold field goes to a finite 
value, £oo, even if v — > oo. This finite value, £oo, is given by the equation, 

£oo (^: arcsin( £oo M) + X ) + ^p(^ 1 ~ (^Mf - cos(M X )) = 0. (5.10) 

Through a numerical calculation, we find that for a choice of x = — vr/M and M = 4, the normalized 
threshold field tends to £oo/er — 0.725 as v tends to infinity. 

The reason why X = — vr/(2M) is critical is that x = — vr/(2M) is the value which determines 
whether the commensurability potential and the impurity potential compete or not at the impurity 
site, namely, the system contains frustration or not. To investigate the frustration, we consider the 
local potential energy Ui at the impurity site, 

UiMxi)) = cos(M^(xi)) - v cos(x + Hxi)). (5.11) 

The first term of the right hand is the commensurability potential and the second term is the 
impurity potential. If the potential energy by the impurity is minimized, that is x + <f>( x i) = 0, 
Ui((f>(xi)) is 

U i(-X) = cos(M X ) - v. (5.12) 
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In the case of —tt/M < x < — 7r /(2Af), the commensurability potential and the impurity potential 
compete and the system contains frustration. In the case of — tt/(2M) < x < 0, however, they do 
not compete and the system contains no frustration. When the system is frustrated the effect of 
the impurity potential to the threshold field is important. However, if the system is not frustrated, 
the impurity potential is irrelevant. 

The reason why the threshold field is finite even if v tends to infinity is that the depinning object 
is not a particle, but a string. In the limit of v — > oo with x = — vr/M which is the optimal case, 
the cusp of the ground state configuration of the phase is on the top of the barrier of the uniform 
potential (Fig. 8). The cusp, however, is just a part of the string, and a finite field is necessary for 



Fig. 8. The classical solution for a choice of M = 4, \ = —ft/M, e = and v — > oo. 



a string to go over the uniform potential barrier. 

§6. The Potential Curve in the Presence of an Impurity 

To clarify the meaning of e c , we consider a potential curve in the presence of an impurity. In a 
homogeneous case the phase is uniform in the ground state, and hence the potential energy density 
(Fig. 3) can be expressed with respect to the uniform phase variable. With an impurity, however, 
the phase is not uniform, and therefore we should consider the potential curve in a functional space. 
In this case the most suitable variable is the phase value, <fn, at the impurity site with x, v and e 
fixed. For each value of <f>i, we can determine the lowest energy configuration; we solve eq. (|3.1[) 
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under the boundary condition of 



and 



(Xi 



1 



M 



1 



e ( — arcsin(eM) -<j )i )+—[^l- ( £ M) 2 - cos(M^ 



(6.1) 



which is similar to eq. (3^). Substituting this optimal configuration into eq. (|2.2| ), we determine its 
energy and obtain the potential curve as the function of (pi which is shown in Fig. 9. While there 
exists an energy barrier below the threshold field, e c , (Fig. 9 (a) and (b)), the barrier disappears 
at e c (Fig. 9 (c)) and CDW depins. 



Fig. 9. The potential curve in the presence of an impurity for a choice of M = 4, \ — —ir/M and v = 1; (a) is in 
the case of e/e c = 0.90, (b) is in the case of e/e c = 0.98 and (c) is in the case of e/e c = 1. 
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§7. Conclusion and Discussion 

We conclude that the existence of an impurity in the model one-dimensional commensurate CDW 
causes the lowering of threshold fields by a finite amount due to the appearance of the local change 
of the phase variable, which is considered as nucleation, near an impurity site. Below the threshold 
field, £t, of the uniform depinning, such local change in the homogeneous case always requires a 
finite excitation energy. However, the nucleation in the presence of an impurity can take place 
without any excitation energy at field, e c , which is smaller than ex- The main reason for the 
lowering of the threshold field is due to the frustration which is caused by the competition between 
the commensurability potential and the impurity potential. 

Our result will be applied to the commensurate CDW with dilute but with macroscopic number 
of impurities where the inverse of the impurity density is smaller than the phase coherence length, 
£. In such a case these will be a distribution of parameter Xi an d then that of the depinning field. 
However the depinning will be triggered by the nucleation at the optimal impurity site, % = —tt/M. 

In actual experiments three dimensionality and quantum fluctuations^ will play important roles 
at low temperatures, which will be studied elsewhere. 

To the best of our knowledge, our investigation is the first microscopic theoretical studies on 
nucleation triggered by the local inhomogeneity in the bulk region. 
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